Self-Consistent Measurement and State Tomography of an Exchange-Only Spin Qubit 
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We report initialization, complete electrical control, and single-shot readout of an exchange-only 
spin qubit. Full control via the exchange interaction is fast, yielding a demonstrated 75 qubit rota- 
tions in under 2 ns. Measurement and state tomography are performed using a maximum-likelihood 
estimator method, allowing decoherence, leakage out of the qubit state space, and measurement 
fidelity to be quantified. The methods developed here are generally applicable to systems with state 
leakage, noisy measurements, and non-orthogonal control axes. 

PACS numbers: 



Nanoelectronics show great promise as a quan- 
tum information platform, in particular as su- 
perconducting qubits [1 5j and spin qubits in 
semiconductors [6 9j. One or two electron spin 
qubits use, respectively, oscillating magnetic [7J 
or electric fields [jSj ilOj . or quasi-static Zeeman 
field gradients |lll414j . to achieve full qubit con- 
trol. Adding a third spin provides exchange- 
driven qubit rotations along two axes, hence full 
control of spin information via electrostatic gat- 
ing only fi", ig EH [13^19]. 

The three-electron exchange-only spin qubit 
has a more complicated level structure than its 
one- and two-electron counterparts [4], I17|, 118] . 
providing, for example, multiple initialization 
states, but also allowing leakage out of the qubit 
state space. Here, we characterize the perfor- 
mance of the three-electron spin qubit by per- 
forming measurement and state tomography [5', 
[20, 21j. Measurement tomography allows accu- 
rate state tomography in the presence of noisy 
measurements and leakage. 

A thrcc-elcctron linear triple quantum dot was 
formed by Ti/Au electrostatic gates patterned on a 
GaAs/AlGaAs heterostructure with the two-dimensional 
electron gas 110 nm below the surface (see Fig. [ija)). 
Left and right plunger voltages, Vi and Vr, controlled 
electron occupation of each dot. All manipulations kept 
a three electron total, with the arrangement, (TV; N^), 
set by the detuning parameter, e = {Vr — V^)/2 — {Vi — 
Vi'^)/2, where e = is defined as the center of 111 (see 
Fig. da)). 

Three electrons have eight possible spin states, four 
with total spin S = 3/2, and four with S = 1/2 [ll22]. 
An external magnetic field splits the eight states into 
four subspaces with spin projection, ms = ±3/2, ±1/2. 



The linear geometry allows two exchange interactions, 
which lower the energy of singlet-like pairs within the 
S = 1/2, ms = ±1/2 subspaces. In particular, tunneling 
between the left and middle dots opens a splitting, Ji{e), 
between the left singlet-like ^ state \Si) = ^(Itlt) - 

litt)) and the left triplet-like state \Ti) = ^(|itt> + 
Itit) - 2 Itti))- Jiis) increases as the detuning is shifted 
towards the 201 charge state. Tunneling between right 
and middle dots similarly opens a splitting Jr(e) between 

\Sr) = ^(itti) - nm and IT,) = ^(itti) + int) - 

2 litt)) which increases as e is shifted towards 102. 

The logical qubit space is chosen to be in the S — 1/2, 
ms = +1/2 subspace [23], where gate voltages control 
the energy spectrum. The logical qubit states, |0) = 

j-, (im) + i;n> - 2 nm and |i) - ^ (im> - i;n», 

are eigenstatcs in the center of 111, with J/(e) — Jris)- 
Two states with 5* = 3/2 couple into the logical sub- 
space through Zeeman field gradients. Longitudinal gra- 
dients couple the qubit space to the S = 3/2, ms = 1/2 
state, \Q) ^ ^ (|m> + int) + nm- The state \Q) 
is a spin symmetric state, being triplet-like for both 
left-middle and middle-right pairs of spins, and is the 
dominant leakage state for this qubit. Leakage into the 
S = 3/2, ms = 3/2 state, |Q-i-) = |ttt), is suppressed 
by a large Zeeman field except at two anticrossings. By 
traversing these anticrossings diabatically — unlike in pre- 
viously triple-dot experiments [HIS] — leakage into \Q+) 
can be made negligible. 

The left singlet-like state, \Si), is prepared by mov- 
ing to ef in 201, and briefly moving near the 201-101 
charge transition border to promote rapid relaxation to 
the ground state. The right singlet-like state, l^,), is 
similarly prepared by moving to in 102 and pulsing 
near the charge border. The excursions to these charge 



FIG. 1: Device, qubit Bloch sphere and spec- 
trum, (a) False color micrograph of lithographically iden- 
tical device with locations of triple dot (smaller red circles) 
and sensor dot (larger red circle). Gate voltages Vi and Vr 
set the charge occupancy of left and right dot as well as the 
detuning, e of the qubit. (b) A Bloch sphere representation 
of the qubit with control axes Ji and Jr indicated, as well as 
two initialization states, \Si) and \Sr)- (c) Energy levels as 
a function of detuning for the lowest energy states |4]. The 
red and blue levels form the logical subspace inside 111, with 
the logical states |0) and |1) denoted at the detuning at which 
they are the eigenstates of the system. Each state has a spin- 
split partner state with opposite spin projection, not shown. 
Values of e for preparation (P) and measurement (M) in 201 
and 102 are indicated. 



borders during initialization are the only departures from 
(5 = 0, where 5 = (K - K°)/2 + {Vi - V°)/2 defines the 
center hne of 111 (see Fig. [2]^ a)). 

Arbitrary qubit states are determined by projection 
onto \Si) or \Sr)- Projection onto \Si) is accomplished 
by moving to the left measurement point, ej^, where \Si) 
can move to 201 while \Ti) and \Q) remain trapped in 
111 During the measurement, an rf excitation is ap- 
plied across the sensor quantum dot. The reflected signal 
is demodulated using homodyne detection [T] and inte- 
grated for tm = 50 /iS, resulting in a signal corresponding 
to either the 201 or 111 charge state. Projection of the 
qubit state onto 15"^) is carried out in a similar way at 
measurement point in 102. 

Effects of exchange interactions, J;(e) and Jr(e), on 
qubit dynamics are modeled by an effective Hamiltonian 

Hj{e) ^ Ji{e)(Ji + Jr{e)crr, (1) 

where ct; = (\/3o-^ -cr2)/4, ci,. ee (-\/3(Ta;-a-^)/4, and 
and tJz are Pauli matrices in the logical basis {|0), |1)}. 
As illustrated in Fig.jljb), J;(e) and Jr{s) drive rotations 
about axes that are 120° apart on the Bloch sphere. In 
what follows, we use the terms J/ and Jr rotations and 
axes in the spirit of this model. 

To demonstrate two- axis control and readout, as well 
as to test the applicability of the simple model, Eq. (1), 
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FIG. 2: Charge stabihty diagram and rotations 
around two axes, (a) Triple dot charge occupancy Nj Nm 
Nr as a function of Vi and Vr in and near the 111 regime; e = 
iVr-V°)/2-{Vi-Vt'>)/2, 5 ^ {K-K°)/2 + {Vi-VP)/2 
The charge occupancy is measured using the change in the 
reflected rf signal, Av^t, incident on the proximal sensor, 
(b) Schematic of a pulse sequence that prepares \Si) in 201, 
and transfers that state to 111, by moving along e at 5 = 0. 
The sequence then waits at for a time rs, and returns to 
201 for measurement. The probability. Pi, of remaining in 
the initial state, 15*;), is plotted as a function of pulse posi- 
tion and wait time. Positive e brings the state closer to 102, 
while negative e brings the state closer to 201. (c) Schematic 
of a pulse sequence, along with a plot of the probability of 
remaining in \Sr} for an excursion to the separation point 
for a time rs. (insets) Model of qubit evolution as a function 
of exchange. No noise has been included. 

we first initialize the system in the \Si) state and sepa- 
rate the electrons into 111 at a detuning for a time rs, 
where the qubit evolves under Hj{e^) [Fig. [2|b)]. The 
qubit is then pulsed to ef^ to measure the projection of 
the evolved state onto \Si), which we determine by mea- 
suring the singlet return probability on the left. Pi, over 
an ensemble of repeated experiments. Pulsing instead to 
allows for a measurement of the projection onto \Sr), 
which when averaged over an ensemble gives the singlet 
return probability on the right, P2. 

Figure |2jb) shows that for states initialized in \Si), 
there is a rapid oscillation of the measured Pi as a func- 
tion of Ts at positive detunings e^, J/(e^) ^ Jr{e^), 
and a roughly constant Pi 1 at negative detunings, 
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FIG. 3: Fast rotation and visibility model, (a) Qubit 
precession as a function of detuning at the separation point, 
for fixed pulse duration, rs — 1.667 ns. This pulse width is 
less than the rise time of both the coax and the pulse gener- 
ator, leading to a significant pulse attenuation at the sample, 
(b) One cycle of the pulse sequence schematic used to gener- 
ate rotation about the Jr axis, (c) Amplitude of oscillation, 
(5 Pi, as a function of dJr/de^, measured from the period of 
oscillations in (a). Theoretical form SPi = ^ ^-'^^(•ij^/'i^^ f ^ 
where a = ts ae/V2 [27| with single parameter, = 450 ^V, 
fit over the first 29 oscillations, to the left of the green point 
(4"' from the left) dot in (a). 



Ji{e^) > Me^). The reverse is true for states prepared 
as IS", ) in Fig.[2jc): i-2 ~ 1 at positive detunings while P2 
exhibits rapid oscillations as a function of rs at negative 
detunings. 

The insets of Fig. [2] show model calcula- 
tions of Pi = I (S'/|e^*-«"-^(^')^s/'*|S'/) |2 and P2 = 
I (Sri e"*-^-^(^'')^s/'* IS'r) p, which agree well with exper- 
iment. These calculations neglect noise in Ji{e^) and 
Jr{£^) as well as fluctuations in local hyperfine fields. 
These contributions are considered in detail below. 

An exchange pulse can generate rapid qubit evolution 
on nanosecond time scales, faster than dynamics induced 
by other sources such as spin-orbit or hyperfine coupling. 
The short-pulse regime thus allows exchange and its noise 
to be examined in isolation from other sources of qubit 
dynamics. Figure 3(a) shows Pi for a short exchange 
pulse, Tg — 1.667 ns, as a function of pulse amplitude, 
Jr{e^), over a range of phase = TsJr{s^)/fi from to 
^ 1587r, corresponding to a 47.4 GHz rotation. 

At large positive e®, where Jr(e^) ^ J/(e^), the noise- 
less model predicts Pi — 5/8 + 3/8 cos(rs Jr(e^)/fi-) for 
initial state \Si) [28 . Experimental data agrees well with 
the 5/8 average [see Fig. 3(a)], but the observed oscil- 
lation amplitude, 5Pi, is notably less than 3/8, with a 
distinct dip where phase varies most rapidly with e^, 
i.e., where dJ^/de^ is largest. The reduced amplitude 
can be understood quantitatively as the result of averag- 
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FIG. 4: Effects of electrical and nuclear noise, (a) 

A schematic of the pulses for a rotation about the Jr axis, 
(b) Free induction decay (FID) as a function of rs. The 
dashed gray curve is a theory model of a rotation around 
the Jr axis (J; = 0) in the presence of dephasing, the solid 
gray curve is a model of Ref. that accounts for the same 
rotation in the presence of dephasing and leakage due solely 
to nuclear magnetic field gradients. The solid black curve 
is fit to a numerical integration of the qubit evolution for a 
rotation around J,., in the presence of quasistatic Gaussian 
distributions of nuclear gradients ABi and ABr, as well as 
quasistatic noise in Jr. Ji was also allowed to be nonzero in 
the fit. Noise is evaluated as quasistatic, rather than incor- 
porating unknown spectral densities. Here, Ji = 20 ± 8 neV, 
Jr = 388±2 neV, and standard deviations as ~ 2.0±0.1 mT, 
aj — 19 ±2 neV, with a turn-on of the exchange modeled as 
an exponential with time constant r — 12.6 ± 0.2 ns. 



ing over exchange noise arising from noise in e, yielding 
dP^ = |e-"'(dJr./deS)=^ ^here a = TsajV2 [17]. The 
period of oscillations in Fig. 3(a) gives a direct measure- 
ment of dJr/de^, leaving a single fit parameter, a^, the 
effective standard deviation of noise in e. Experiment 
and theory are in excellent agreement [Fig. [3](c)]. The fit 
value, (Te — 450 /iV, is only nominal, as it includes ef- 
fects of finite coax bandwidth, making it larger than the 
actual e noise in the system. 

In Fig |4j the free induction decay measured at longer 
separation times, rs, and at fixed detuning, e^, reveals 
the combined effects of exchange noise, which causes de- 
phasing, discussed above, and Zeeman field gradients, 
which cause both dephasing and leakage out of the qubit 
space. 

Quasi-static longitudinal (effective) field differences 
between dots, ABi = Bf - and ABr = - B^, 
drive coherent evolution between |0) and |1). Gradients 
due to hyperfine fields appear static on the time scale of a 
single sequence of measurements, but execute a thermal 
random walk over an ensemble of measurements. In ad- 
dition, Zeeman differences g^,^ABi and g^^AB^ compa- 
rable in magnitude to Ji{£) or Jr{s) will drive evolution 
into IQ), the leakage state. Here /^b is the Bohr mag- 
neton and g ~ —0.4 is the electron g factor. Averaging 
over the entire nuclear ensemble during repeated mea- 
surements results in a damped oscillation towards the 
triplet outcome as the qubit dephases and leaks into \Q). 
By examining in detail the rg dependence of Pi at a fixed 
in Fig. |4j and comparing it with theoretical models 
for low frequency exchange and Overhauser [3] noise, we 
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conclude that nuclear fluctuations are the predominant 
source of noise in this system, with a standard deviation 
of 2.0 mT. 

Low-frequency hyperfine noise 30, 31 can be compen- 
sated using dynamical decoupling |32l |33j . Unlike the sit- 
uation in double quantum dots, however, a single-pulse 
echo cannot undo the effects of two hyperfine field gra- 
dients in the three-dot systeni[19j. Nevertheless, a single 
TT-pulse can undo a portion of the dephasing due both to 
nuclei and low-frequency exchange noise. 

Single-pulse partial echo is demonstrated by preparing 
\Si) in 201, separating to 111 where the state rotates for 
a time ti around J^, followed by a 7r-pulse around J;, fol- 
lowed by further rotation around for a time T2. The 
sequence is illustrated in Fig. [sja) . Pi shows robust os- 
cillations as a function of both sum and difference of the 
dephasing times ti and T2 , similar to a Ramsey measure- 
ment with a refocusing pulse in the middle. The decay 
envelope in ti -|- T2 at ri = T2 in Fig. [sjjb) gives a lower 
bound on the coherence time, T2 ^ 100 ns, while the 
decay envelope in ti — T2 shown in Fig. [Sjjc) gives a de- 
phasing time, T2 ^ 25 ns. The dephasing time in ri — T2 
is consistent with FID times [Fig. |4|^b)], while decay as 
a function of ri -|- T2 is extended by a factor of ~ 4 for 
the echo condition ti — T2. This modest enhancement 
is consistent with decoupling a portion of the noise from 
the environment. A model of classical, slowly fluctuating 
hyperfine field gradients [Fig.jsjjd)], yields a value for the 
standard deviation of ABi and AB^. of 3.4 mT, and in- 
dicates these gradients to be the dominant noise source 
for this pulse sequence. 

Tomographic Characterization of the System 

Qubit performance is commonly assessed by state and 
process tomography. These techniques require a well- 
characterized set of measurements that give enough infor- 
mation to reconstruct all the matrix elements of the den- 
sity matrix of the system. It is important to recognize, 
however, that measurements themselves are subject to 
noise, relaxation, and systematic errors, and may not be 
well described by idealized projective measurements. To 
accommodate both state preparation and measurement 
errors we implement a self-content approach that com- 
bines measurement data and models of system dynamics 
within a maximum likelihood estimation routine. Self- 
consistent tomography along similar lines has been car- 
ried out recently for superconducting qubits in Ref. [S]. 

Consider single-shot measurement in the singlet-triplet 
basis, for the moment ignoring leakage into state \Q). 
The singlet measurement fidelity Fs is the probability 
that the measurement correctly registers the singlet out- 
come when measuring a system prepared in the singlet 
state; similarly. Fx is the probability that system pre- 
pared in the triplet state yields the triplet outcome for 
the measurement [2] . Given that the measurement is de- 
fined to be along the \Si)-\Ti) or \Sr)-\Tr) axis of the 
Bloch sphere, these two numbers completely character- 
ize an imperfect two-outcome measurement on a qubit. 
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FIG. 5: Dynamical decoupling, (a) A schematic for a 
three pulse sequence demonstrating an echo, where a pre- 
pared IS';) processes around Jr for a time ri, performs a n- 
pulse around Ji, and then precesses again around Jr for a 
time T2 before being measured in 201. (b) Probability Pi of 
measuring \Si) for the pulse sequence depicted in (a), for the 
case ri = T2. The solid black curve is a model with noise 
from Jr, ABi, and ABr during the dephasing times ri and 
T2, and noise due to Ji, ABi, and ABr during the 7r-pulse 
using parameters extracted from FID data in Fig. |4|b). (c) 
The results of a three pulse echo sequence, illustrated in (a), 
that maps the probability of remaining \Si} as a function of 
the total dephasing time (ri + T2), and the difference in time 
between the first free induction decay, n, and the second, T2. 
(d) A model plot for (c) that averages over thermal distri- 
butions of nuclear gradients. See Supplementary Information 
for details of the calculation. 



On the other hand, measurements along other directions 
require rotations which are themselves imperfect as well. 
This, then, requires a general description of a noisy mea- 
surement that includes errors in measurement direction 
as well as reduced fidelities. Such a description is pro- 
vided by the formalism of Positive Operator- Valued Mea- 
sure (POVM) elements. As all our measurements have 
two outcomes, the POVM describing each measurement 
basis choice i is given by a single positive Hermitian ma- 
trix Ei associated with the "singlet" outcome. (The cor- 
responding "triplet" outcome is associated with the ma- 
trix (I—Ei).) The eigenvectors of Ei determine the axis of 
the Bloch sphere along which the measurement is made. 
The eigenvalues are bounded between zero and one, with 
the larger eigenvalue of Ei equal to Fs while the smaller 
eigenvalue is 1 — Ft- Using the POVM formalism, the 
probability that a system described by a density matrix 
p will yield the singlet outcome when the measurement 
of basis i is performed is 



P,{p) = Ty[E,p]. 



(2) 



The POVM formalism can be applied to the three- 
state system of qubit states plus leakage state \Q), in 
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which case the Ei are 3x3 Hermitian matrices. As our 
measurements are insensitive to coherence between the 
leakage and qubit states, which in any case is expected 
to be small, we restrict Ei to be incoherent with the \Q) 
space (i.e., each Ei has support on the reduced qubit sub- 
space together with a population in \Q)). We note that 
the \Q) state will, with an idealized spin-to-charge mea- 
surement, always yield the triplet outcome. As such, the 
I Q) -population of Ei quantifies the error-induced proba- 
bility that a system prepared in the \Q) state will instead 
yield the singlet outcome. 

Measurement Tomography 

Measurement of the four matrices, Ei, i=l-4, is re- 
quired for state tomography, additionally yielding the 
population of the leakage state. Each of the four Ei 
has five parameters associated with the qubit state 
and the leakage population, for a total of 20 unknown 
quantities |35j . Measurement tomography therefore 
requires measurement outcome statistics on five well- 
characterized input states, pj, using each of the four gen- 
eralized measurements to yield 20 independent observed 
probabilities, Pij. We then solve P^- = TT[EiPj] for Ei 
self-consistently, subject to the constraints on the eigen- 
values of Ei to lie between and 1. Details are given in 
the Methods section and the Supplementary Information. 

Using the reconstructed Ei, we extract the measure- 
ment bases and fidelities from the eigenvalues and eigen- 
vectors as discussed above. The singlet outcome fidelity 
is indicated by the length of the Ei arrow in Figs. |6je,f). 
The smaller eigenvalue in the qubit space relates to the 
fidelity of the triplet outcome, Fxi = 1 — Xi2, while the 
eigenvalue in the leakage subspace relates to the proba- 
bility that \Q) will have a triplet outcome, Fgi = 1 — A^s. 
Finally, the measurement visibility in the qubit subspace 
is found as Vi — Fsi+FTi — 1 = Xii — Xi2- The average sin- 
glet fidelity over all four generalized measurements was 
found to be 69%, while the average triplet fidelity was 
80%, giving an overall average measurement fidelity in 
the qubit subspace of 75% with an average measurement 
visibility of 49%. 

State Tomography 

We can now perform state tomography on arbitrary 
states of our system using our set of tomographically 
characterized generalized measurements despite the fact 
that these measurements are inherently noisy. As a 
demonstration, we generated sets of unknown states by 
performing a simple rotation around Ji in the presence of 
dephasing for multiple fresh input states. In Fig. [6j two 
separate input states, pi (red) and p2 (blue), are prepared 
and then pulsed to a negative detuning where Ji ^ J^, 
followed by a generalized measurement (the markers in 
Fig.[6];a-d)). Using our descriptions of Ei , we are able to 
reconstruct the state at a set of time intervals during the 
evolution by solving equation ^ again for p (the mark- 
ers in Fig. |6][e,f)). The theory curves overlaid on top of 
the data in Fig. Iwa-h) are fits to the Liouville-von Neu- 
mann equation, accounting for the finite bandwidth of 



the coax and a theoretical model of the exchange profile, 
in the presence of the nuclear noise determined from the 
calibration procedure. 

As expected, pulsing towards a negative detuning 
yields an outcome that depends on the input state. For 
the states prepared as pi = \Si){Si\, the red markers and 
curves in Fig. |6]^e,f), sitting at a position of large Ji only 
imparts a trivial phase. At this detuning, \Si) is split 
off energetically from \Q), suppressing leakage out of the 
qubit space. This is observed on the Bloch sphere as a 
collection of points near the idealized \Si) state. 

The state that is prepared as p2 — \Sr){Sr\, the blue 
markers and curves in Fig. [6|e,f), has a very different 
response. Since it is an eigenstate of Jr, at this detun- 
ing it is a superposition of the ground and excited qubit 
states. As a result, it precesses around the J; axis in the 
presence of dephasing, which causes the state to spiral 
inwards towards the rotation axis. 

At this detuning a fraction of p2 is in the excited qubit 
state, which is energetically close to the \Q) state, al- 
lowing the Overhauser gradients to rotate that fraction 
out of the qubit space. This leakage into the \Q) state 
occurs on a 10 ns timescale in both the data and the 
theory, which increases the decay of the Bloch vector r 
towards the center of the Bloch sphere in Figs |6f^e,f,h) 
as the probability exits the qubit subspace. The leakage 
is seen clearly in the rise of the \Q) population, Pg, in 
Fig. [BJg). The qubit state vector length in Fig.[6]jh) acts 
as a measure of qubit coherence and population, and it 
decays with a 16 ns, which is consistent with pre- 

vious measurements. 

In conclusion, we have demonstrated initialization, 
two-axis electrical control, and self-consistent state re- 
construction of an exchange-only spin qubit. The ex- 
change interaction allows extremely fast qubit operation. 
The method of tomographic calibration we developed can 
be applied directly to other qubit systems to quantify 
measurement errors. Future work will include investi- 
gating regimes where J; and Jr are simultaneously much 
larger than the nuclear gradient Zeeman energy, which 
would suppress leakage into \Q), and structures compris- 
ing of six dots that implement a fast two-qubit gate jl5) . 
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Methods 

Measurement tomography requires a choice of five in- 
put states, pj, which span the qubit subspace. The ini- 
tialization states 15*;) and \Sr) provide pi and p2 respec- 
tively. We create two additional states p3 and p^ by 
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FIG. 6: Measurement and State Tomography. State tomography of the qubit during evolution around Ji. Measurement 
outcome probabihties for four measurement bases for initial states pi (red) and p2 (blue), pulsed towards near the 201-111 
charge transition, producing rotation mostly around Ji. The solid curves in (a-d) are a fit to the model (see Supplementary 
Information). Solid curves in (g) and (h) are generated from the model in (a-d). (a) Schematic and measurements of Pi(pi) 
(red) and Pi(p2) (blue); (b) P-iipi) (red) and P2(p2) (blue); (c)P3(pi) (red) and P3(p2) (blue); (d) P4(pi) (red) and P4(p2)(blue). 
(e-f) Views of Bloch sphere with measurement axes. Graphical representation of the qubit portion of E\ (red), E2 (blue), 
(green) and £4 (black), (g) Population of the leakage state as a function of separation time, (h) The length of the Bloch vector 
as a function of separation time data (blue triangles), model (blue curve), and gaussian fit [width 16.4 ± 0.9 ns] (dashed blue 
curve). 



rotating \ Si) around Jr at e| and \Sr) around Ji at ef re- 
spectively. These four input states span the qubit space. 
The rotated states are subject to rotation errors, dephas- 
ing, and leakage, which we need to characterize using a 
phenomenological model for the dynamics and fit the pa- 
rameters of this model using experimental data. 

To facilitate this, we use a series of states for and 
P4. The noisy evolution of p3 and is then modeled 
with a generalization of Hj{e) to the larger manifold of 



|0), and IQ), including the effects of AB/ and AS,, 
(see Supplementary Information for details). During the 
calibration of the EiS, we compare the model of this evo- 
lution to the series of states P3(r3) created by initializing 
\Si) and rotating around at e| for a set of times T3 be- 
fore measuring. The series of states P4(t4) was produced 
in a similar fashion by preparing \Sr) and rotating for 
a set of times r4 about Ji at £4. These series of states 
contain enough information to determine the strength of 
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the nuclear dephasing and the exchange axes at ef and 
£4. 

The final input state, ps, is chosen to be a completely 
mixed state with no coherences remaining and a signifi- 
cant weight in the leakage state. This choice allows for 
accurate measurements of the parameter in each 

Ei. The ensemble of ps was prepared by performing re- 
peated pulses to dephase around J; and Jr over a distri- 
bution of ABi and ABr- 

With the observed statistics Pij for five known input 



states j using four measurements i, we determine Ei by 
fitting the calibration probabilities from all of our input 
states to our model of the noisy evolution to produce a 
Maximum Likelihood Estimate (MLE) for E1-E4 as well 
as the nuclear noise and exchange during the calibra- 
tion |36]. For the data in Fig.|6j the standard deviations 
of ABi and AB^ were ~ 2.5 mT, which is consistent with 
the earlier estimations extracted from the FID and echo 
data. 
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Supplementary Information for Self-Consistent Measurement and State Tomography of an Exchange-Only 

Spin Qubit 

This supporting document describes further details of the fabrication, state readout, noise modehng, and measure- 
ment tomography techniques. The measurement tomography section details the pulse sequences and fitting routines 
used to extract the POVM elements, as well as the effects of finite bandwidth limitations on the state reconstruction. 



The three-electron system was confined in a lateral triple quantum dot formed in the two-dimensional electron gas 
(2DEG) at the GaAs/Alo.sGao.yAs interface 110 nm below the surface of the heterostructure. The GaAs/Alo.aGao.yAs 
heterostructure was grown on a solid-source Varian Gen II molecular beam epitaxy (MBE) system equipped with an 
arsenic valved cracker source to provided As^ for the growth. The heterostructure was grown on a semi-insulating 
(100) GaAs substrate with a growth rate of 1 /im/hr. The 2DEG is formed by a Si modulation doping ((5-doping) of 
~4 X 10^^ m~^ (40 nm away from the 2DEG interface). Hall effect measurement done at 20K gives a 2DEG density 
of ^2.6 X 10^^ m~^ and a mobility of ^-^43 m^/V s. 

High bandwidth coaxial lines were attached to the left, middle, and right plunger gates of the triple quantum 
dot, and a radio-frequency (rf) reflectometry circuit was connected to a neighboring quantum dot for fast state 
readout [3 [5] . The experiment was performed in a dilution refrigerator equipped with a cryogenic amplifier (noise 
temperature T/v ~ 3 K), with an electron temperature of ~ 120 mK. An in-plane external magnetic field of 300 mT 
was applied along the dot connection axis [see Fig 1(a) in main text]. 



A uniform normalization procedure was used for all data in Figs. 2-5 to convert the measured reflectometry signals 
into output probabilities. For a given set of pulse parameters (e^, Tg, T1+T2, T1-T2, etc.), the qubit was measured 
using four preparation and measurement routines. In the first two routines, the state 15;) was prepared in 201 then 
measured either in 201 {\Si) projection, yielding Pi) or in 102 (l^^) projection, yielding ^2)- In the other two routines, 
the state \Sr) was prepared in 102 then measured either in 201 (\Si) projection) or in 102 {\Sr) projection). Each 
measurement consisted of sitting at the measurement point — ef^ (e^) for \Si) {\Sr)) readout — and integrating the 
demodulated rf signal reflected from the impedance transforming circuit [1] attached to the rf-sensor quantum dot for 
tm — 50 fis to yield v^[. 

This process was then repeated while stepping one of the pulse parameters (e^, rg, T1+T2, T1-T2, etc.). Each sequence 
was then repeated 2^^ or 2^^ times to obtain measurement statistics. The resulting data was then histogrammed, 
following the procedure in Ref. and fit to a function of the form. 



A. 



Device 



B. 



Measurement and Normalization 



1. Normalizations Based on Single- Shot Outcomes 




0.00' 
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Figure SI: Single-shot histogram. A histogram of outcomes for Fig. 3(a) in the main text. Red solid curve is a fit to 
equation (|STJ. 
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(SI) 



where ni^v^-f) is the fraction of histogram events with outcomes Wrf for a measurement in 201, v^f^ is the reflected 
voltage corresponding to double occupancy in the left dot, is the reflected voltage corresponding to single charge 
occupancy in the all three dots, Av^i = v^l^ — v^f^, Pi is the fraction of 201 outcomes in the data set, Ti is the 
relaxation time at ej^ , tm is the measurement time, and a is the standard deviation of the histogram peaks due to 
noise in the rf equipment and shot noise intrinsic to the rf sensor dot. For measurements in the right dot, rir has an 
identical form, with all 201 notations replaced with 102 and Pi replaced with P2. 

The extracted parameters ^2^^^ and are then used to normalize the return probabilities on the left side. 



(S2) 



where {vrf{s^,Ts)) is the average voltage over all repetitions of the measurement sequence for a specific and rg. Pj 
is normalized similarly, with v'^f^ replaced by v^^^. 



Equation ( S2 1 converts v^f into a probability, but it does not account for relaxation during the measurement time 
Tm, where a 111 state relaxes to a 201 state for 15*;) projections or a 102 state for 15*^.) projections. Relaxation during 
the measurement was accounted for using a two step process. The histogram shape is only weakly dependent on the 
precise value of Ti, but failing to allow for relaxation of the 111 charge state would underestimate the separation 
between histogram peaks. Ti decay was therefore included in equation (SI I to determine the peak positions v'^f^ 



and v^l^ accurately, but the Ti fit parameter is not itself an accurate measurement of the relaxation time in the 
data. In order to more accurately correct for relaxation, we project a state prepared as \Sr) in 102 against \Si) in 
201, and record the probability as Pf°', at the beginning of each sequence and compare it with the theoretical value 
|(S'r|S'i)P = 0.25; we confirm the theoretical value by measuring relaxation as a function of measurement time tm, as 
described in Sec. B.2. Traces were then corrected as 



(S3) 



A measurement of a state prepared as \Si) and measured in 102 is similarly used to correct for Ti decay in P2. 



2. Ti Relaxation During Measurement 



The relaxation time was extracted from a measurements of the dependence of the uncorrected singlet probabilities 
Pi" and P2 as functions of the integration time tm for \Si) measured in 102 and 15^) measured in 201. The observed 
dependence was well described by exponential relaxation of triplets integrated over the measurement time, 



1 - P" = — / dt{l - A)e~'/^' 
Tm Jo 



n (1 - A)Ti , 
po ^ 1 - ^ ^(1 



Tm 



(S4) 



where A is the fraction of singlets present at tm = 0, and corresponds to |(S'r|5'z)P in the absence of any evolution 
when pulsing the prepared state to the measurement position. The expression for P^ is identical, with Ti referring 
to the measurement relaxation on the right rather than the left. In the data shown in Fig. |S2[ extrapolation to zero 
measurement time yields A = 0.24 ± 0.02 for both Pj* and P2 . This allows us to correct for relaxation by normalizing 
the data with the theoretical singlet probability. 
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Figure S2: Measurement relaxation and confirmation of 120° axis separation (a) Schematic for preparing \Si) and 
projecting it onto \Sr}- (b) Experiment (triangles) and numerics (solid curve) for the probability of measuring a singlet in 
102 if the state was prepared as jS'i) in 201, as a function of measurement time tm- The numerical results are uncorrected for 
relaxation during the measurement. Numerical result is a fit to the function (S4| with fit parameters Ti — 137 ± 9 ps and 
\{Sr\Si)\' = 0.24 ± 0.01 at tm = 0. (c) Schematic for preparing \Sr} and projecting it onto \Si). Experiment (triangles) and 
numerics (solid curve) for the probability of measuring a singlet in 201 if the state was prepared as \Sr) in 102, as a function of 
measurement time tm, uncorrected for relaxation during the measurement. Numerical result is a fit to the function |S4| with 
110 ±8 Ats and \{Sr\Si)\^ = 0.24 ±0.02 at tm =0. 



fit parameters Ti 



C. Figure 4b Theory Curves 



1. Pure Electrical Dephastng 



We model exchange rotation in the presence of low-frequency detuning noise as: 



-Pi(rs) g ^ 1 ~ cos[rs J(ts)] exp 



(S5) 



where T2 is a characteristic time for dephasing due to electrical noise. Rise-time effects, due for instance to bandwidth 
limits of the coaxial cable, are modeled as an exponential rise in the exchange, 



J{t) = Jo 



1 — exp 



t 



(S6) 



The dashed gray curve in Fig. 4b of the main text used the parameters Jq = 391 neV, Tj* = 25 ns, and tr = 13 ns. 



2. Dephasing due to Nuclei - Ladd Curve 



Following Ref. [3], Eqs. 9, 14-16, we model the effect of nuclear fluctuations on the qubit as 



A(t-s) = ^ + ^^-/i(ts) + cxp 



/2(rs) 



(S7) 
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(S8a) 
(S8b) 
(S8c) 

(SSd) 



Here, we take the standard deviation of the nuclear fluctuations to be the same in each dot, which simplifies the 
expressions in Ref. '3'. The solid gray curve in Fig. 4b of the main text uses fit parameters Jo = 386 ± 2 neV, 
CAB = 1-9 ± 0.2 mT, and tr = 13 ns. 



3. Numerical Model: Electrical and Nuclear dephastng 

In order to incorporate both electrical and magnetic sources of noise, we used a numerical model for the time 
evolution of an initial 15;) in the presence of both exchange interactions Ji and J^, as well as longitudinal field gradients 
ABi and ABr- The finite bandwidth of the coax and function generator are accounted for with an exponential turn-on 
as described in Sec. C.l. We take Jr to be Gaussian distributed, appropriate for small amplitude fluctuations in 
over a range where varies approximately linearly with e^, that is, 6 J « (AJ^/ Ae^)Se^ . In the region of detuning 
where dJ^/de^ 3> dJj/de^, only fluctuations in Jr were taken into account. Since the left exchange was decreasing 
while the right exchange was increasing, they were approximated as: 

Ji = e-*/""^ (S9) 
= J^{l - e'*/^") + SJ (SIO) 

The effects of the slowly fluctuating nuclear bath were incorporated by taking an ensemble average over Gaussian 
distributions of nuclear gradients between the left and middle (ABi) and middle and right (ABr) dots. Limiting 
this model to detuning regions away from the \Q+)-\Si) and \Q+)-\Sr) anticrossings, transverse components of the 
hyperfine field can be safely neglected, leaving only gradients between longitudinal components. This model assumes 
a Gaussian distribution of classical nuclear gradients, with no back-action on the nuclei from the qubit. 

Explicitly, Pi was evaluated numerically using a uniform step size, 

Piirs) = I '^^^['l^Y^'' {Si\e-^^^^-'^\Si) %-(As?+AB^)/(2.^)-(.J)V(2.J) (Sii) 



(271)3/2^^^ 



X ^ {ABjyidx ABjji^f^) {Jmax Jmin) 



where To, is the time spent during the rotation, gb is the standard deviation of the Gaussian distribution of classical 
values that each nuclear gradient could achieve, aj is the standard deviation of the Gaussian distribution of Jr values, 
Nstep is the number of discrete values sampled for each Gaussian, ABmax and ABmin are the limits of ABi and ABr 
values sampled, Jmax and Jmin are the limits of 5J values sampled, tr is the turn-on time for the exchange. Hi is 
the Hamiltonian at the dephasing position. The Hamiltonian consisted of two parts, the model laid out in Ref |4j, 
and a nuclear Hamiltonian; 'H(e^) — HjiJi, Jr) + ^'Hb{ABi, ABr), where 7 ~ g^iB = — 25.4neV/mT. 
We can write the exchange Hamiltonian in the basis of (|0)-|1)-|(5)) as: 

(-liJl+Jr) ^{Jl~Jr) 0\ 

^{Jl^Jr) -\{Jl+Jr) (S13) 

\ 00/ 
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Here, the zero energy state has been shifted relative to Hj{s) in the main text to make the energy of the \Q) state 
zero at zero detuning. This brings our notation into agreement with Ref. [4]. 
The longitudinal nuclear terms in this basis are: 



( U^Br 
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hABr) 
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(S14) 



where ABi = {Bf — and ABr — (^1 ~ -^1) differences in local magnetic field along the z direction. Terms 

that only contribute a global phase in this basis have been dropped for clarity. 

A fit to this model yields: = 21±8 neV, J° = 388±2 neV, as = 2.0 ±0.1 mT, aj = 19±2 neV, tr = 12.6 ±0.2 
ns. The distributions were each sampled evenly 12 times (Ngtep) each for a total of 12'^ — 1728 samples between 3aB 
and — 3ctb for the nuclei and between 3aj and —Saj for 6 J. This gives the solid black curve in Fig. 4b. 

D. Figure 5d Echo with Hyperflne Dephasing and Leakage, without Electrical Noise 



The partial echo (Fig. 5 in the main text) was analyzed using a model similar to the one used in Figure 4c. The 
timescales involved in the echo are much longer than the rise time tr, so the phenomenological exponential turn on 
of the exchange is removed for simplicity. In addition, since the dephasing was dominated by nuclei in Figure 4 the 
noise on is omitted. 

The model includes the evolution of an initial \Si) under the action of the three exchange pulse sequence in the limit 
of instantaneous rise times in the qubit environment of 111. The pulses are evaluated in a piecewise-static manner, 
ignoring the weak adiabatic effects associated with pulsing from one detuning position to the other. Explicitly, Pi 
was evaluated numerically using a uniform step size as: 



E 



dABi dAB,. 



-(ABf+AS^)/(2a|) 



-(AB?+AB^)/(2<T^) 



(S15) 



(S16) 



where ti(2) is the time before (after) the tt pulse, gb is the standard deviation of the Gaussian distribution of nuclear 
gradients, Nstep is the number of discrete values sampled for each gradient, ABmin and ABmax are the limits of 
ABi and ABr values sampled, 'H(e^) is the Hamiltonian at the dephasing position, includ ing n uclei, and 'H{£^) 
is the Hamiltonian for the tt pulse. The Hamiltonians were of the same form as equations (S13) and (S14), with 
l-i{e^) = nj{0,Jr) + -fl-LB{ABi,ABr), ^(e^) = nj{Ji,0) + jHsiABi, ABr). The values used in the model were 
7 — gfiB = — 25.4neV/mT, as = 1.7mT, Jr = 276 ncV, J; = 824ncV. The magnetic field gradients were each 
sampled uniformly between 3ctb and — 3ctb, with Nstep = 40. 



E. Measurement Tomography for the Exchange Only Qubit 

1. Measurements and Measurement Operators 



In order to determine the populations (|0)(0|, |1)(1|) and coherences (Rc[|0)(l|], lm[|0)(l|]) in the qubit subspace 
(four unknowns), as well the population of the leakage state (|Q)(Q| = 1 — |0)(0| — |1)(1|), we need to perform four 
measurements. The measurement probabilities can be expressed in the following fashion: 



P,{p) = Tr[E,p] 
P2{p) = Tr[£;2p] 
Psip) = Tr[E3p] 
Pi{p) = TriE^p] 



(S17a) 
(S17b) 
(S17c) 
(S17d) 



where p denotes an unknown input state and Ei is a measurement operator that describes the fidelity of a singlet 
outcome for a measurement in the i*"^ basis If we have a set of five (or more) known input states, pj, one can use 
those states to measure the Ei, E2, E^, and by solving the set of equations (S17). Once the Ei are determined, 
we can reconstruct any unknown state p from the four probabilities. Pi, P2, P3, and P4. 
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Our measurement tomography approach uses known input states to characterize the measurement operators Ei. 
Of the five required input states, three are relatively easy to prepare. These are the two initialization states, \Si) and 
\Sr), and the completely mixed state. Their density matrices in the |0)-|1)-|(5) basis can be written as: 



pi = \Si){Si\ - U# i (S18) 




P2 = \Sr){Sr\ = i (S19) 





1 

3 

.0 



P5 = \mixed){mixed\ = I \ I (S20) 



We create ps by pulsing to regions of large J; and repeatedly, allowing the state to dephase around both rotation 
axes as well as the nuclear gradients. The preparation of the dephased state is confirmed by comparing measurements 
of the dephased state that was initially prepared as \Si) with the state initially prepared as \Sr)- The only way that 
these two outcomes will be identical is if they are both completely dephased in the qubit space. The length of the 
sequence is many times larger than Tj^^^, which when combined with the pulses yields a completely mixed state 
with \Q) as wellf^. This state allows for the characterization of the \Q) fidelity in each measurements, which is not 
necessarily identical to the qubit triplet-like {\Ti) and [Tr)) fidelity. 

Those three states are entirely real by construction. It is more challenging to prepare high fidelity states with 
lm[|0)(l|] ^ 0, which are needed to fully characterize the system. These superposition states allow us to characterize 
the complex quanties of our measurement operators, and in turn allow us to measure any superposition of states 
in the qubit subspace. Since we do not have access to initialization states with lm[|0)(l|] ^ 0, we need to create 
these states through evolutions under control pulses. These pulses themselves contain noise. The extent to which we 
correctly account for the dephasing and leakage that occurs during the preparation of these states determines our 
ability characterize the two rotated input states, — Rr\Si){Si\Rl and p4 — Ri\Sr){Sr\Rj , and therefore E3 and 
E4. We incorporate the details of the evolution, including noise, into the estimation of our input states using MLE 
techniques, which results in a more accurate state reconstruction. 

The final algorithm we use simultaneously estimates the evolution of pa and p4 as well as the POVM elements 
El — Ei by using a Levenberg-Marquardt least-squares method to minimize the difference between the output of the 
model and the measured probabilities associated with pi — P5. This gives us our Maximum Likelihood Estimate for 

El — i?4. 

2. Creation and Determination of the Known Input States 

One way to minimize the error in estimating ps and p4 is to spread ps and p4 out into a set of input states that 
evolved under a common axis of rotation for a range of rotation times T3 and T4. In other words, pa becomes a set of 
states P3(t3) evolving at and p4 becomes a set of states P4(r4) evolving at s"^. With these sets of rotation times, we 
can accurately estimate axis of the rotations and the frequency of the rotation at the fixed detunings e'^ and and 
therefore each state in the series P3(t3) and p4(r4). In this manner we can prepare a set of reasonably high fidelity 
input states that are superpositions of our measurement basis states with complex coefficients. 



Explicitly, we take the effective Hamiltonian for this system, equations (S13l and (S14|, and evolve the initial 



density matrix in time. We then average that evolution over a distribution of AB's, calculate the probabilities using 



equations (S17) and compare it to the data. The fitting procedure, if started from good initial guesses (which we 
generate iteratively through incremental fits, and incremental relaxation of constraints), can generate the MLE for 
the Hamiltonian parameters and every element of Ei, E2, E3, and E4. 

The fit is performed by calculating the theoretical Hamiltonian of the system, using the same techniques that were 
employed for the FID and the partial echo. Assuming time independent Hamiltonians, P3(t3) and P4(t4) can be 
written as 

J 27r(T^ 



P4(r4) = / d^|i^e-«(-')-p2e^«(-')-e-(^^?+^^^)/(2'^^), (S22) 
ZTTag 
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Figure S3: Measurement tomography pulse sequences. Schematics for pulse sequences that create the five input states, 
(a) Input state pi. (b) Input state p2. (c) Input state ps prepared from an initial 15*;). (d) Input state pa. (e) Input state p4. 
(f) Input state ps prepared from an initial \Sr)- 



where ^(e^) = Hj(J,(e'), Jr(e')) + 7^i3(AS,, AS,), ^(e^) ^ nj{Me^), Me^)) + jneiABi, ABr), and <jb is 
the standard deviation of the Gaussian distribution of nuclear gradients. We then generate the probabihties for 



each measurement using Eqs. (S17). The Levenberg-Marquardt least-squares algorithm subsequently optimizes the 
parameters of 'H{e^), 7^(£^), and (Tb to minimize the difference between measured and calculated probabilities. The 
estimates that the fitting function produce come with error bars, which may be useful indications of the reliability of 
the MLE output. 

There is a further improvement that we can make to this scheme. Ps,{ts) and Pa{ts) are formed by pulsing to a 
region of high exchange, and then pulsing to the settle point, while pi and p2 are formed by pulsing from 201 and 
102 respectively. This can increase the lowpass effects that the settle point is trying to mitigate. To improve the 
situation, and to have pi, p2, and ps have equal weight with P3(ts) and P4(ts), we can evolve pi and p2 under the 
same Hamiltonians as pa{ts) and P3(ts) respectively, where pi is approximately an eigenstate of the Hamiltonian that 
evolved P4(ts), and p2 is approximately an eigenstate of the Hamiltonian that evolved P3(ts). 

A few further improvements were made to the calibration routine to mitigate the low-pass effects of the coaxial 
lines; 

• A voltage overshoot was added to the first 833 ps of the pulse to for P3(t3) and p4(r4). This makes the pulse 
shapes closer to the ideal square pulse, which is easier to evaluate in our fitting routine. 

• The calibration routine records the evolution of P3(t3) and ^4(14) for T3,r4 > 4 ns, which further reduces the 
effects of transients at the beginning of the pulses. 

• The settle point employed directly before measurement is a tradeoff between decoupling the rotation pulses from 
the measurement pulses, which improves with settle time tw, and the dephasing brought on by nuclei, which 
gets worse with increasing r-yv. The settle time needed increases with the amplitude of the pulse directly before 
the settle point. In order to keep tw to a minimum, pi is placed at e** for ti = 5 ns, the evolution point for 
/04(t4), where it is approximately an eigenstate. This prevents the need to pulse all the way from 201 to the 
settle point, which is a significantly larger amplitude pulse. p2 is similarly placed at for T2 = 5 ns to reduce 
the amplitude necessary to bring it to the settle point as well. The time spent at e'^ and is then incorporated 
into the fitting routine as well, making Pi(ti ~ 5ns) and P2{t2 — 5 ns). 

• — 5ns), P2{t2 — 5ns), ps are then repeated to give them the same weight in the fitting routine as P3(t3) 
and ^4(74). 

As the last refinement, we allow for small preparation infidelities, which are included as incoherent mixtures of 
\Ti) + \Q) and |T'r) + |Q) into Pi(ti — 5ns) and P2{t2 = 5ns) respectively. This adds two more free parameters. 

There are 27 unknowns in this problem: five unique quantities for each of the four measurement operators, two 
unknown exchange terms at and two at e^, one unknown for the standard deviation of the nuclear gradients, and 
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two unknowns for the preparation infidelities. These 27 parameters were optimized over the entire data set, resulting 
in the best estimate of the POVM elements and the Hamiltonians at the calibration positions. 

Once the system is fully characterized, we invert equation (S17l, this time using known i?i's, to solve for the 
unknown p. This solution is then used to create the MLE output using the techniques described in Ref. [6]. 



3. Normalization of Single-Shot Data in Figure 6 



Unlike the charge sensor normalization procedure used for the data in Figs. 2-5, which was covered in Sec. [B] 
the data in Figure 6 was not normalized by extracting v^^^ and and normalizing the average voltage. Instead, 
we used a single-shot threshold voltage, w^, which was chosen to separate 201 and 111 outcomes, as was done in 
Ref. [2]. This allowed us to make a clean comparison to previous single-shot measurements, including the definitions 
of measurement fidelity. The downside of this method is that it reduces the visibility of the oscillations, which is part 
of the reason for the reduced amplitude in figure 6. Pi is then defined as the fraction of outcomes whose v-^f is on the 
201 (for Pi and P3) or 102 (for P2 and F4) side of v'^^. 





Figure S4: Measurements of pj. The calibration curves and their fits, (a) p\ (b) p2 (c) ps (d) pi (e) ps 



Figure [S4] shows the example data for the calibration routines which determine E1-E4. The markers are the 
measurements of the known input states pi-p5, while the solid curves are the outputs of the MLE routine that 
extracts the Hamiltonian parameters and the POVM elements. The calibration for pi, p2, and p^ are repeated 
to provide an equal weighting with p^ and p4 in the MLE routine. The Hamiltonian parameters are: for ps^T^): 
Jiie^) = 33 ± 8 neV, Me^) = 380 ± 4 neV, for P4(t4): Ji(e'') = 530 ± 8 neV, Me'^) = 4 ± 10 neV, as = 2.8 ± 0.1 
mT. These parameters define four measurement axes, which when inverted yield the data in Figure 6. 



F. Theory Curves in All Panels of Figure 6 



In the presence of sharp (instantaneous) pulses, the theory curve for the \Sr) initial state stays close to the plane 
that cuts though \Sr) and the origin, perpendicular to J;. This is because the pulses move rapidly through the center 
of 111, where Ji and Jr are both active. If that region is traversed slowly, the net rotation pulls the state towards 
jiSi), with the end result being that subsequent rotations around Ji trace out a larger circle on the Bloch sphere as 
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it approaches a great circle. To correct for this increased circumference, the fuU bandwidth hmited pulse shape is 
needed. 

The theory curves in Fig. 6 of the main paper were generated by solving the Liouville-von Neumann equation, 
ihp = ['^(e), p], where 'H(e) is a function of the time-dependent pulse detuning, as well as the nuclear gradients, ABi 
and ABj.. Ji and are modeled following Ref. with the modification that the tunnel couplings are themselves a 
Gaussian function of the detuning. This modification was necessary to explain the data in Figs. 3 and 6. 
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(S23) 



(S24) 



Here, the tunnel coupling width Wt = Sq/O.SS mV, the charge transition detuning 
a = 40/ieV/mV, and the tunnel coupling t = 22 /leV. 



18 mV, the lever arm 
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To model the system accurately, the ordinary differential equation (ODE) was solved up to the end of the rotation 
at for each value of rs- Those solutions were then used as the initial values in a second solution, where the system 
was pulsed to the settle point for a time tw before pulsing the the measurement point. The pulses were modeled as 
the low passed output of the intended piecewise functions, 

r£° + (£i-e")[l-exp(-tlo)] iitQ<t<ti 
s^t) = \{^° + i^' " - exp(-%)]) exp(-*^) (S25) 

where is the measurement point, is a resting point at large exchange prior to the rotation, and is the rotation 
point . The time spent at the resting point starts at to — —2.5 ns, while the time at the rotation point starts at 
ti = ns. The set of pulses used in the second ODE were 

Seit)^ {eArs)expi-'-^) if rg < t < rg + rw, (S26) 

which accounted for the approach to the settle point. After this, the POVM elements take effect, describing the rest 
of the measurement. 

The time evolution is calculated for each configuration of ABi and ABr, and averaged over a Gaussian weighting 
for each gradient, the standard deviation of which was extracted in the POVM characterization routines described 
above. No exchange noise was included in this model. The measurement probabilities associated with the averaged 



solution of the Liouville-von Neumann equation were then extracted using equations (S17|. The population of the 



\Q) state was extracted as Pq = The Bloch sphere vector length is extracted from the qubit subspace of the 

density matrix as \f\ — ^ {p*xY + {PyY + (Pz)'^: where 



_ (mo\ io)(i|\ . . 

p [\im \i){i\) ^^^^> 

is the density matrix of the qubit subspace. Here, p* = Tr[p*(T3;] = 2Re[|0)(l|], p* — Tr[p*CTy] = — 2Im[|0)(l|], and 
p:-Tr[pV;] = |0)(0|-|l)(l|. 
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